Anisotropic imbibition on surfaces patterned with polygonal posts 
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We present and interpret lattice Boltzmann simulations of thick films spreading on surfaces patterned with 
polygonal posts. We show that the mechanism of pinning and depinning differs with the direction of advance, 
and demonstrate that this leads to anisotropic spreading within a certain range of material contact angles. 



I. INTRODUCTION 



A drop placed on a partially wetting substrate will make a finite angle with the surface, given by Young's equation 0, 
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where 7s v> 7sl and 7 are the solid-vapour, solid-liquid and liquid-vapour surface tensions. Young's equation assumes that the 
surface is smooth, and that the contact line is able to move freely to allow the drop to globally minimise its free energy. 

Due to recent advances in microlithography it is now possible to pattern surfaces with regular arrays of micron-scale posts, 
leading to deviations from Young's equation on a macroscopic level. On a superhydrophilic, or superwetting surface, the fluid 
can be drawn into the spaces between the posts, such that the drop forms a film with thickness equal to the height of the posts 
a phenomenon that is termed imbibition. 

Imbibition is thermodynamically feasible if the thick film has a lower free energy than the dry surface. The free energy change 
per unit width, S T, when the film advances a distance 5x can be estimated by averaging over the surface features 



SJ 7 = [(jsl ~ 7sv) (r - 4>) + 7 (1 - </>)] Sx 



(2) 



where r is the ratio of the surface area to its vertical projection, </> is the fraction of the surface covered by posts, and we assume 
posts of constant cross- section. Eliminating the surface tensions in Eqn. ^ using Eqn. ([I]), the condition 5F < becomes Q 



COS#y > COS#i 



(3) 



This inequality relies on the same assumption as Young's equation ([T]): namely that the contact line can move freely over the 
substrate, sampling the average properties of the roughness. This assumption holds well on some surfaces, for example, in the 
longitudinal direction on a grooved surface, but not in other cases, for example perpendicular to such grooves, where free energy 
barriers due to contact line pinning can halt the motion of the interface 0. 






(a) (b) 

FIG. 1: (a) Illustration of the Gibbs' Criterion. At a sharp corner a wetting contact line remains pinned for all angles between By and By + ip 
(between the full lines), (b) Diagram showing the geometry and post dimensions for the simulations described in Sec. [in] 
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Contact line pinning occurs when an interface moving across a surface meets a convex corner. The criterion for pinning, 
proposed by Gibbs | 6 ] and demonstrated experimentally by Oliver et al. (7), is that the contact angle can take a range of values 
spanning the dihedral angle of the corner, as shown in Fig. |l(a)| Over this range, the contact angle with respect to the dry plane 
is too low for the contact line to advance, and that with the wet plane is too high for the contact line to recede. Pinning on surface 
features can lead to the threshold angle for imbibition being substantially lower than predicition ([3]), as was demonstrated by 
Courbin et al (SI EJ - Furthermore it has been shown that anisotropic surface features lead to anisotropic spreading |[T0lfTT1l . In 
these proceedings, we use the lattice Boltzmann method to study imbibition through an array of posts with uniform polygonal 
cross-section, building on previous work fT2l . We show that the mechanism of contact line pinning differs with direction, and 
explain how this leads to anisotropic spreading. 



II. SIMULATION APPROACH 



We model the sytem as a diffuse-interface, two-phase fluid in contact with a solid substrate. The thermodynamic state of the 
fluid is described by an order parameter p(r), corresponding to the density of the fluid at each point r. The equilibrium properties 
are modelled by a Landau free energy functional over the spatial domain of the fluid V, and its boundary with solid surfaces dV, 

* = Iff K - 2 £ r wd - ^ 2 ) - 1} - + H vp| 2 ) dv - [[ fjL sP ds . (4) 

The first term in the integrand of ^ is the bulk free energy density, where v = (p — p c )/p c and p c ,p c , and j3r w are constants. It 
allows two equilibrium bulk phases, liquid and gas, with v = ivT^w- The second term is a Lagrange multiplier constraining 
the total mass of the fluid. The third term is a free energy cost associated with density gradients. This allows for a finite- width, or 
diffuse, interface to arise between the bulk phases, with surface tension 7 = | p c Y // 2^p c (/3r w ) 3 and width x = \ Pc v^/(/^ r wPc)- 
The boundary integral takes the form proposed by Cahn [13). Minimising the free energy leads to a Neumann condition on the 
density 

<9±p = -/is/ft • (5) 
The wetting potential p s related to the Oy of the substrate by ifTTl 



p s = 2 /3t w \J 2p c ftsign (f — Oy) y cos § (l — cos ^) , a = arccos (sin 2 Oy) • (6) 
The hydrodynamics of the fluid is described by the continuity and the Navier-Stokes equations 

d t p + d a (pu a ) = , (7) 
dt(pu a ) + dp(pu a up) = -dpP a p + dp (prj {dpu a + d a up} + p\S af3 d 7 u 7 ) , (8) 

where u is the local velocity, P is the pressure tensor derived from the free energy functional ^ and r\ and A are the shear and 
bulk kinematic viscosities respectively. A free energy lattice Boltzmann algorithm is used to numerically solve Eqns. ( |7|8| ) 1 14 - 
fT6l . At the substrate we impose the boundary condition ^ E7][T8]|, and a condition of no-slip lfT9ti2l1 . 

We choose k = 0.01, p c = 0.125, p c = 3.5, rw = 0.3 and ^ = 1.0, giving an interfacial thickness x = 0-9, surface tension 
7 = 0.029 and a density ratio of 3.42. The viscosity ratio is t^l/^g = 7.5. 



III. IDENTIFYING THE PINNING MECHANISMS 



We consider a rectangular array of posts on a flat substrate. The cross-section of each post is uniform, and is an equilateral 
triangle, oriented to point along a primary axis of the array taken to be the x-direction (see Fig. 1(b) ). In our simulations we hold 
the array spacing d and post side-length b at 40 and 20 lattice units respectively and vary the post height h. We find that rescaling 
the system such that d = 60 does not change the threshold angles of spreading significantly 1 23 ]. The posts and substrate are 
taken to have the same Young angle 0y- 

We consider the advance of a straight contact line, which is parallel to the ?/-axis. We exploit periodic boundary conditions and 
use a simulation box of length d y along y. We further halve the computational burden by taking x = as a plane of reflectional 
symmetry, and we compare the dynamics for triangles pointing away from, or towards, the origin. To simulate imbibition fed by 
a mother drop resting on the surface would require a very large simulation box, and be prohibitively costly in terms of computer 
time. Since we are only interested in the details of flow amoungst the posts, we instead feed imbibition from a 'virtual reservior', 
a small region ~ 6 lattice points wide spanning the centre of the box where v is fixed to y/ (3t w at each time step of the simulation. 
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FIG. 2: Snapshots of (a) the connected contact line and (b) the disconnected contact line mechanisms of pinning. In the middle image of (b) a 
gap appears between the interface and the face of the post. This is because the liquid (blue) surface represents the density p c , and close to the 
concave corner the wetting potential increases the density above p c 

In this way, liquid is introduced while there is outwards flow, but once the interface is fully pinned, no new liquid enters the 
system. 

6y is then decreased quasistatically, and we record the value at which depinning and spreading to the next post occurs. For the 
geometry we describe, Eqn. ([3]), which describes imbibition with no pinning, gives an upper bound of the threshold angle of 
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A. Pinning of a connected interface 

Snapshots showing one pinning mechanism, for a film advancing in the direction of the points of the triangles, are shown in 
Fig. 2(a) The film is of height h up to the leading triangle, and then descends with increasing x to meet the substrate at the 
Young angle. There are two ways in which the contact line can move forward. Firstly, it could make a shallower angle at the 
substrate, but this would increase the free energy away from the minimum characterised by Eqn. ([I]). Secondly, the top of the 
film could move forwards, but this would create liquid-gas interface and hence also have a free energy cost. 

Depinning will occur when Oy is sufficiently small that the contact line on the base reaches the next post. This depinning 
pathway, which we shall term the connected contact line mechanism, was elucidated by Courbin et al[ 8 , 9 ], who showed that 
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Numerical results for the variation of the depinning angle with h/b are shown in Fig. [3] as indigo circles, and Eqns. d9l [To]) are 



FIG. 3: The threshold angle for depinning along +x (indigo circles) and — x (mauve squares) as a function of h/b. Depinning in these 
directions occurs according to the connected and disconnected contact line mechanisms respectively. The predictions from Eqns. (|9) in red 
(dotted), {To) in blue (dashed), and |T2] > in green (full) are plotted for comparison. 



plotted as red and blue curves respectively. Comparing the simulation data to the blue curve, we see that the simulation values 
are significantly higher than those predicted, expecially for lower values of h/b. 

To resolve the discrepency we note that Eqn. (10) assumes a flat interface. A positive Laplace pressure Ap will instead 
produce a convex curvature, enabling the interface to extend further across the substrate. Neglecting curvature in the y direction, 
we model the interface in the xz plane as a circular arc with radius of curvature R = j/Ap given by Laplace's law. The contact 
angle with the substrate will then be modified to 



tan(#i-/3) = f , 

where f3 is the angle of bulge, given by \Jh? + s 2 = 2R sin /?. The depinning threshold is thus given by 



0i = arctan [7] + arcsin 
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We expect the dominant contribution to the Laplace pressure to result from confinement in the z direction. Therefore, we shall 
assume Ap oc h~ x . Writing R = Ah and s = d — Bb, a least squares fit of the data to Eqn. ( 12 ), with respect to A and B, was 
performed. The optimisation found B = 0.822, barely different from the value ^ « 0.866 used in Eqn. ([9]), and A = 7.09. 
Eqn. ( [T2| ) is plotted with these coefficients as the green curve in Fig. [3] and the fit is very reasonable. 



B. Pinning of a disconnected interface 

We now present simulation results with the posts pointing towards the origin, and identify a second mechanism for (de)pinning, 
shown in Fig. |2(b)| Now the advancing front is disconnected, and is pinned at the vertical edges of the posts. The base of the 



film is pulled forward by the hydrophilic substrate, but there is a free energy cost associated with the growth of the interface as 
it spreads out from the gap. As the Young angle is quasistatically decreased, the contact line creeps onto the blunt faces of the 
posts, near to the base substrate, but remains pinned to the post edges at higher z, where the angle made between the interface 
and the blunt faces remains less than 0y. When 0y becomes sufficiently small, the depinned parts of the contact lines from 
neighbouring gaps meet each other midway. Once connected, the interface readily wets up the posts and out across the substrate. 

We shall refer to this as the disconnected contact line pinning mechanism. The threshold for depinning is plotted in Fig. [3] as 
mauve squares. The dependence on h/b is different to that for motion along +x. When h/b is low, the depinning angle closely 
follows the upper bound given by Eqn. ([9]), indicating that pinning by the posts is weak in this regime. For larger values of the 
ratio h/b, 6j levels off to - 51°. 
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FIG. 4: Spreading on a square array of triangles 






(b) 



(c) 



FIG. 5: Hexagonal posts show the connected contact line mechanism along ±x and the disconnected contact line mechanism along ±y. 

IV. IMBIBITION THROUGH POLYGONAL POSTS 



We now present simulation results for films spreading through arrays with various lattice symmetries and post geometries. 
We use arrays which are several posts wide in both the x and y directions, such that the film is not connected over periodic 
boundaries. We again use a virtual reservoir, this time located at a small location at the centre of the array, but we hold Oy 
constant over time. We discern how both the arrangement, and the geometries, of the posts affect the dynamics of the interfaces, 



and the final film shapes. These can be interpreted in terms of the pinning mechanisms identified in Sec. [111 

For our simulations we use d 
inhibit spreading in the +x and 

times in the evolution of the film. The posts are shown in brown, the wetted substrate in blue, and the unwetted substrate in 
white. 



40, b = 20, h = 30 and Oy = 55°. According to Fig.[3j these parameters should allow and 
x directions respectively. In Figs. [4] [5] and [6] we show plan views of the substrate at various 



A. A square array of triangular posts 



We first consider the system studied in Sec. Ill extended in the y direction. The shape of the film at intermitant times is shown 
in Fig. [4] Advance of the film is possible in the +x and ±y directions, via the connected contact line mechanism, but the film is 
barred from advancing in the — x direction, where the disconnected contact line mechanism, which has a lower threshold angle, 
is relevant. Thus the surface acts as a microfludic diode. Such unidirectional behaviour is made possible by the triangular shape 
of the posts. 



B. A square array of hexagonal posts 

Having considered exclusively triangular posts thus far in the chapter, we now turn our attention to posts whose cross-sections 



are regular hexagons. We find that the two depinning mechanisms discerned for triangles, in Sec. |III| may also be applied to 
hexagons, but that their directional distribution of occurrence is different. We consider hexagonal posts in a square array, oriented 
so that the corners point along ±x. Along these two directions, as might be expected, the (de)pinning behaviour follows the 
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FIG. 6: Spreading on a hexagonal array of traingles 



connected contact line mechanism. Conversely, the faces point along ±y, and it is the disconnected contact line mechanism 
which determines the (de)pinning in these directions. Fig. [5] shows the spreading of a film on a square array of hexagons. Since 
advance of the liquid is permitted along ±x but barred along ±y, a stripe of fluid is formed. 



C. A hexagonal array of triangular posts 



We now simulate a hexagonal lattice of posts with spacing d = 40, and the triangles aligned with lattice directions, as shown in 
Fig. [6] We start with a circular film with diameter spanning several posts (Fig.[6ja)). As spreading begins, the film quickly facets 
into a hexagon, by aligning its sides with posts in the immediate vicinity (Fig. [(Jb)). Spreading continues, via the connected 
contact line mechanism, along the directions of the three corners of the posts, but the interface is pinned, by the disconnected 
contact line mechanism, along the faces. As a result, the facets along the corner directions shrink as they advance (Fig.[6jc)). 



V. DISCUSSION 



We have performed Lattice Boltzmann simulations of imbibition on hydrophilic substrates patterned with posts, whose cross- 
sections are regular polygons. Our motivation was to identify pinning mechanisms on the posts and show how these lead to 
anistropic spreading behaviour on the surface. 

We began by considering the advance of a long planar front along a row of triangular posts. This enabled us to take advantage 
of periodic boundaries in the simulations, reducing computational expense, and to isolate particular pinning behaviours. The 
simulations showed that the critical value of By at which the interface advances differs between directions relative to the triangles. 
Hence there is a range of By in which spreading is unidirectional, with the exact range and direction of the anisotropy depending 
on the relative dimensions of the substrate. The cause is differing depinning routes: one where the contact line along the base 
substrate is connected, and one where it is disconnected, punctuated by the blunt edges of the posts. 

We showed that a square lattice of triangular posts inhibits spreading in one direction, while if hexagonal posts are used, the 
spreading is bidrectional, with films elongating. Finally we investigated spreading amongst a hexagonal lattice of triangular 
posts. The three-fold rotational symmetry of this geometry leads to the formation of a triangular film. 

In future work it would be of interest to consider how the spreading is affected if the post cross section changes with height, 
and how electrowetting might be used to locally control the contact angle, and hence the spreading characteristics (22). 
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